Shear-induced fragmentation of Laponite suspensions 
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Simultaneous rheological and velocity profile measurements are performed in a smooth Couette 
geometry on Laponite suspensions seeded with glass microspheres and undergoing the shear-induced 
solid-to-fluid (or yielding) transition. Under these slippery boundary conditions, a rich temporal 
behaviour is uncovered, in which shear localization is observed at short times, that rapidly gives 
way to a highly heterogeneous flow characterized by intermittent switching from plug-like flow to 
linear velocity profiles. Such a temporal behaviour is linked to the fragmentation of the initially 
solid sample into blocks separated by fluidized regions. These solid pieces get progressively eroded 
over time scales ranging from a few minutes to several hours depending on the applied shear rate 
7. The steady-state is characterized by a homogeneous flow with almost negligible wall slip. The 
characteristic time scale for erosion is shown to diverge below some critical shear rate 7* and to 
scale as (7 — 7*)~" with n ~ 2 above 7*. A tentative model for erosion is discussed together with 
open questions raised by the present results. 

PACS numbers: 83.60.La, 83.50.Rp, 83.60.Pq, 83.50.Lh 



INTRODUCTION 

Interest into the solid-to-fluid transition displayed by 
"soft glassy materials" (SGM) has risen tremendously 
in the past two decades not only because of its indus- 
trial importance but also due to its inherent theoretical 
issues and experimental challenges. On the theoretical 
side, the fact that this transition, also referred to as 
the "jamming-unjamming" transition, can be triggered 
either by increasing the temperature, by lowering the 
system concentration, or by applying a strong enough 
external load has led to propose a universal "jamming 
phase diagram" Although this picture remains ap- 

pealing, it is clear that a complete understanding and 
modelling of glassy-like phenomena involved in materi- 
als with microstructures as diverse as submicron hard 
spheres, micronsized grains, short-range attractive par- 
ticles, or highly charged platelets, are still out of reach 
Most recent experiments performed to investigate 
the structure and the dynamics of SGM at rest have re- 
lied on single or multiple light scattering techniques and 
have focused on ageing properties M or on the presence 
of dynamical heterogeneities 0, @, 13- Together, ageing 
and dynamical heterogeneities call for a spatiotemporal 
modelling of SGM that is still lacking. 

Even more difficult is the task of elucidating the flow 
mechanisms in these out-of-eqnilibrium systems when 
submitted to some external shear. Rheology has long ap- 
peared as the only tool to probe the mechanical response 



'Present address: Brandeis University, Physics Department, 415 
South Street, Waltham, MA 02453, USA 



of SGM and the shear-induced solid-to-fluid transition, 
also called the "yielding transition." In particular, the 
problem of defining and measuring correctly the yield 
stress ay, i.e. the critical shear stress below which the 
material behaves as a solid and above which it flows like 
a liquid, has focused a lot of attention The diffi- 
culties raised by the yield stress most often result from 
thixotropic features due to the competition between age- 
ing and shear- induced "rejuvenation" of the sample [9|. 
Thus estimating precisely ay involves applying or mea- 
suring very small shear rates over waiting times that can 
reach several hours. 



Yield stress fluids and shear localization 

Some authors have proposed to distinguish between a 
"static" yield stress ay, defined for the material at rest, 
that would be the analogue of a static friction coefficient, 
and a smaller "dynamic" yield stress ay, measured on the 
flowing material by decreasing the imposed stress, simi- 
lar to a dynamic friction coefficient [lO|. Such a discrep- 
ancy linked to hysteresis effects also raises the question 
of whether the yielding transition in SGM is continuous, 
i.e. the shear rate 7 continuously increases from zero 
when the shear stress a is increased above ay, or rather 
discontinuous, i.e. 7 jumps to some finite critical value 
7c as soon as ct > Cy. 

Continuous models for the flow curve cr vs 7 include the 
Bingham and the Herschel-Bulkley models, which have 
been recognized to hold for many SGM [HI . [T3 | . On the 
other hand, a discontinuous yielding transition questions 
the basic assumption of homogeneous flow that underlies 
any standard rheological measurement [l^, ■ Indeed, 
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if a shear rate 7 < 7c is apphed, the flow is expected to 
display shear localization, i.e. the coexistence of a solid- 
like region where the local shear rate 710c vanishes and a 
fluid- like region where 710c = 7c- 

The recent development of original experimental tools, 
such as active or passive microrlicology 14, [H, light 
scattering under shear, and time-resolved local velocime- 
try, that can go beyond standard rheology, has allowed 
one to address some of the above mentioned issues with 
renewed interest. In particular, using direct visualiza- 
tion, nuclear magnetic resonance (NMR), particle imag- 
ing, or ultrasonic velocimetry, shear localization was ob- 
served on a large variety of SGM, ranging from emulsions 
0, [13, or colloidal suspensions [11 [l^, M, [H, [13 
to wet granular materials 24 , 2^ . Still, in most of these 



experiments performed in Couette geometry, it remained 
unclear whether shear localization could be attributed 
to a truly non-zero 7c or to the stress inhomogeneity 
inherent to the concentric cylinder geometry. It is not 
until very recently that shear localization has been also 
evidenced in the cone-and-plate geometry where stress 
heterogeneity is minimized, thus establishing flrmly the 
relevance of discontinuous models 1261 . 



Structure and rheology of Laponite suspensions 

Within the last decade, among the huge variety of 
SGM, Laponite, a synthetic clay of the hectorite type 
made of heavily charged disc-shaped particles of diame- 
ter 25-30 nm, thickness 1 nm, and density 2.5 g.cm""^, 
has emerged as a good, albeit complicated, candidate 
to explore the above issues. When dispersed into wa- 
ter at a few weight percents, typically 0.6-4 wt. %, 
and depending on the ionic strength, Laponite sus- 
pensions evolve from a low-viscosity liquid to various 
solid-like "arrested" states. The sol-gel transition in 
Laponite has been investigated mostly using light scatter- 
ing 27 , H, 2^, 3^, 31 1 . The ageing properties of Laponite 
have triggered lots of effort and debate about the exact 
nature of these nonergodic states, in particular, about 
whether the material is actuall y a "gel" or a (Wigner) 



"glass" [30|, 132|, l33|, l3J, |35|, l36|, l37 



. Moreover, due 
to their out-of-equilibrium nature, Laponite suspensions 
were used to test possible violations of the fluctuation- 
dissipation theorem (FDT) [3^, \^ . Several recent mi- 
crorheology experiments tend to prove that the effective 
temperature cannot be distinguished from the bath tem- 
perature, so that the FDT remains valid [U, [4^, but 
such results are still under debate (isj . 

The nonlinear rheology of Laponite suspensions has 
also been intensively studied with emphasis on thixotropy 
and ageing or rejuvenation under shear 27, 0, 0, 0, 
47l . H El, [13, EH • A few experiments have hinted to 



the presence of shear localization in the vicinity of the 
yield stress in Laponite samples [3, 5^, 5^, which pro- 



vided support for discontinuous models of the yielding 
transition [l^, [2^. However, to the best of our knowl- 
edge, the local velocity field of Laponite suspensions was 
measured only in a wide-gap (2 cm) Couette cell with 
a temporal resolution of about 25 s per velocity profile 
using NMR velocimetry [H^l and in a plate-plate geome- 
try of gap 7 mm through dynamic light scattering (DLS) 
[53| . The two corresponding papers only showed a couple 
of velocity profiles, the former raising the issue of strong 
geometry-induced stress heterogeneity and the latter re- 
porting both wall slip and slow temporal oscillations of 
the velocity with the important limitation that DLS only 
provides point-like velocity measurements and necessi- 
tates to mechanically shift the cell in order to scan the 
whole gap. Therefore a systematic timc-resolved study of 
velocity profiles in sheared Laponite suspensions would 
certainly provide new insights into the flow mechanisms 
involved in yielding. 



Summary of our previous and present work 

In a previous work [H^l, we have explored for the 
first time the infiuence of boundary conditions on yield- 
ing in Laponite samples by simultaneous rheology and 
time-resolved ultrasonic velocimetry. For the purpose of 
measuring local velocity using ultrasound, a significant 
0.3 wt. % amount of hollow glass spheres of mean di- 
ameter 6 fiui was added to our 3 wt. % Laponite sus- 
pensions. By comparing two experiments performed un- 
der similar imposed shear rates in smooth and rough 
(sand-blasted) Couette geometries on timescales of order 
5000 s, we unveiled a dramatic effect of surface rough- 
ness on the fiow mechanism during yielding. Indeed, 
while the scenario observed with rough walls was con- 
sistent with a discontinuous transition characterized by 
shear localization in agreement with pre vious observa- 
tions on Laponite suspensions lol . 52, 5^ and on other 
SGM [il, [13, [12, [H, [iil, slippery boundary conditions 
in the smooth cell led to a thoroughly different picture. 
When wall slip was allowed, the sample was reported to 
break up into macroscopic solid pieces that are slowly 
eroded by the surrounding fluidized material up to the 
point where the whole sample has become fluid. 

The aim of the present paper is to provide a full data 
set, more discussion, and a toy model of this original 
yielding scenario triggered by slippery boundary condi- 
tions in Laponite suspensions seeded with microspheres. 
Our paper is organized as follows. Section describes the 
materials and methods used to study the yielding transi- 
tion in Laponite. Results obtained in a smooth Couette 
cell on a 3 wt. % Laponite sample seeded with 1 wt. % 
microspheres are presented in Section for a large range 
of applied shear rates. Finally, Section provides a dis- 
cussion of the results, in particular about the possible 
influence of the microsphere concentration, and an at- 
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tempt to model the experimental observations. 



EXPERIMENTAL 
Sample preparation 

Laponitc powder (Rockwood, grade RD) is used as re- 
ceived and dispersed at 3 wt. % in ultrapure water with 
no added salt. The pH of the solution is approximately 
10. According to the most recent litterature [i^jlsll, such 
a suspension (3 wt. % Laponite without salt) is supposed 
to fall into the Wigner glass region of the phase diagram. 

In the present study, the Laponite suspensions are im- 
mediately seeded with hollow glass spheres of mean di- 
ameter 6 /xm (Sphericel, Potters) at a weight fraction 0.3 
or 1 wt. %. As explained below in Section , these micro- 
spheres act as contrast agents for ultrasonic velocimetry. 
Most of current protocols used for Laponite preparation 
recommend to filter the samples with a 0.4-0.8 /im mesh 
size after vigorous stirring for 15-30 min. Such filtra- 
tion is believed to break up Laponite clusters that may 
form due to incomplete dissolution and may lead to erro- 
neous interpretations of the sample structure as "gel-like" 
from light scattering data [11]. Here, however, due to the 
presence of the microspheres, subsequent filtration of the 
samples is not possible. 

Within 30 min of magnetic stirring, the dispersion be- 
comes homogeneous and very viscous but remains fluid. 
Due to the presence of the microspheres, the solution is 
slightly turbid, which will allow for direct visual inspec- 
tion of the sheared samples. The sample is then left to 
rest at room temperature. After about two hours, the 
sample viscosity has increased by several orders of mag- 
nitude and the dispersion has clearly built up a yield 
stress and solid-like properties since it is able to sustain 
its own weight. We let ageing proceed for at least two 
days, i.e. for waiting times tu; > 2 10^ s. For such large 
iu), we expect the ageing dynamics to have slowed down 
so much that the influence of ageing can be neglected 
on the timescales of typically one hour involved in the 
present experiments. In any case, as discussed below, 
the samples are pre-sheared before any experiments in 
order to erase most of the sample history through shear 
rejuvenation. 



Rheological protocol 

Linear and nonlinear rheological properties are mea- 
sured using a stress-controlled rheometer (Bohlin C-VOR 
150, Malvern Instruments) in Mooney-Couette (concen- 
tric cylinder) geometry. Our cell is made out of Plexiglas 
without further treatment. The bob (inner cylinder) of 
radius i?i = 24 mm is rotating and will be called the "ro- 
tor" in the following. It is terminated by a cone with an 



angle of 2.3°. The stator, i.e. the fixed cup (outer cylin- 
der), has a radius R2 = 25 mm. The gap width of this 
Mooney-Couette geometry is thus e = R2 — Ri = 1 mm. 
The curvature of our cell leads to a stress decrease of 
about 8 % from the rotor to the stator. The height of the 
rotor is = 30 mm. The standard deviation of height 
profiles of the cylindrical walls measured using atomic 
force microscopy is typically 15 nm, which will be re- 
ferred to as "smooth." The whole cell is immersed in a 
water tank of volume of about 1 L connected to a water 
bath whose temperature is kept constant and equal to 
25±0.1°C (sec Fig. [J). 

Before any measurement, the sample is pre-sheared for 
1 min at +1500 s"^ and for 1 min at -1500 s ^. Such a 
protocol erases most of the sample history through shear 
rejuvenation and ensures that the strain of the sample 
accumulated during loading into the cell has no influence 
4J| . We then proceed with a standard oscillatory test at 
1 Hz in the linear regime for 2 min. This test allows us 
to make sure that the values of the viscoelastic moduli 
G' ~ 500 Pa and G" ~ 25 Pa no longer change signifi- 
cantly at the beginning of the actual experiment. We also 
checked that this procedure leads to reproducible results 
over a few hours. Finally, at time t ~ 0, the experiment 
proceeds either with the standard (linear or nonlinear) 
rheological tests shown in Section or with combined rhe- 
ological and velocimetry measurements (see Section ). 
In the following, a and 7 denote the shear stress and the 
shear rate indicated by the rheometer. The rotor velocity 
vq and 7 are linked by: 



i?i(i?i + R2) 



Rl 



Rl 



■76, 



(1) 



where the geometrical factor i?i(i?i -I- R2)/[R\ + R\) ac- 
counts for the cell curvature [H^ . In the presence of wall 
slip or heterogeneous flows, 7, the so-called "engineering" 
or "global" shear rate may strongly differ from the local 
shear rate which will be noted 710c or 7(r), where r is the 
radial distance from the rotor. 



Ultrasonic velocimetry and optical imaging 

Our setup for combined rheology and local velocime- 
try is sketched in Fig. [TJ The sample velocity field is 
measured using ultrasonic speckle velocimetry (USV) at 
about 15 mm from the cell bottom. USV is a technique 
that allows one to access velocity profiles in Couette ge- 
ometry with a spatial resolution of 40 /xm and a tempo- 
ral resolution of 0.02-2 s depending on the shear rate. 
It relies on the analysis of successive ultrasonic speckle 
signals that result from the interferences of the backscat- 
tered echoes of successive incident pulses of central fre- 
quency 36 MHz generated by a high-frequency piezo- 
polymcr transducer (Panametrics PI50-2) connected to 
a broadband pulser-receiver (Panametrics 5900PR with 
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200 MHz bandwidth). The speckle signals are sent to a 
high-speed digitizer (Acqiris DP235 with 500 MHz sam- 
pling frequency) and stored on a PC for post-processing 
using a cross-correlation algorithm that yields the local 
displacement from one pulse to another as a function of 
the radial position r across the gap. One velocity profile 
is then obtained by averaging over typically 1000 suc- 
cessive cross-correlations. Full details about the USV 
technique may be found in Ref. [s^l- 

The scattered signal is provided either by the material 
microstructure itself or by seeding the fluid with "con- 
trast agents" with the constraint to remain in the sin- 
gle scattering regime. In the case of Laponite samples, 
which are transparent to ultrasound, we used the micro- 
spheres described above to produce ultrasonic echoes in 
a controlled way. The sound speed in our samples was 
independently measured to be c = 1495 ± 5 m.s""'^ at 
25°C. 

The microspheres also lead to a slight turbidity which 
allows for a direct visualization of the sheared samples. 
A simple CCD camera (Cohu 4192) is set in front of the 
water tank and the Couette cell is lit from behind. Im- 
ages of the samples are recorded at a frame rate of 5 fps 
and later synchronized with the velocity profile measure- 
ments as explained in our previous work (Hij . In this 
paper, we shall only show images typical of the various 
stages of the experiment. 



RESULTS 
Conventional rheometry 

Standard rheological measurements performed on a 
3 wt. % Laponite suspension seeded with 1 wt. % mi- 
crospheres are shown in Fig. [51 The viscoelastic moduli 
displayed in Fig. and measured in the linear regime 
for a shear stress amplitude of 5 Pa (which corresponds 
to a strain of at most 0.13 % over the accessible range of 
frequencies) are typical of an arrested (gel or glass) state. 
The elastic modulus G" is always much larger than the 
loss modulus G" and only decays from 500 Pa to 400 Pa 
over frequencies / = 0.07-8 Hz. G" also slowly decreases 
from about 30 Pa to 15 Pa over this range of frequen- 
cies and shows no sign of downturn at low frequencies, a 
feature that is usually attributed to the presence of very 
slow relaxation modes in SGM [l2. 14|. 

Figure EJb) shows the evolution of the viscoelastic 
moduli in the nonlinear regime for shear stress ampli- 
tudes above 15 Pa at a fixed frequency of 1 Hz. These 
data clearly point to a yield stress of 47 Pa, when de- 
fined as the shear stress amplitude ay for which G' {ay) = 
G"{ay). This yield stress corresponds to a strain of 25 %. 



All these measurements are fully consistent with previ- 
ous data on Laponite suspensions [32 
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the addition of microspheres influences the ionic strength 
of the suspension. The conductivity of a 1 % wt. suspen- 
sion of our glass microspheres in water was measured to 
be about 100 //S.cm^^, which is roughly equivalent to a 
salt (NaCl) concentration of 0.8 10~^ mol.L"^. Increas- 
ing the microsphere concentration induces a noticeable 
stiffening of our samples: G' increases by about 20 % 
when the microsphere concentration is increased from 
0.3 wt. % to 1 wt. % probably due to the change in ionic 
strength [33 |. 

As explained in the introduction, another way to probe 
the yielding transition using rheology is to perform non- 
linear measurements where a constant shear stress or 
shear rate is imposed. Since the flow curve is almost 
fiat for small shear rates, it is usually more suitable to 
work under imposed shear rate. Figure [U^c) shows the 
fiow curve obtained through a shear rate sweep from 0.3 
to 2500 s~^ within 200 s followed by the correspond- 
ing downward sweep at the same rate. The most ob- 
vious feature of this flow curve is the large hysteresis 
between upward and downward sweeps. This is typical 
of thixotropic materials whose microstructure is slowly 
modified by shear resulting in a time-dependent appar- 
ent viscosity 77 = cr/7. The very same kind of hysteresis 
cycle was observed recently for similar shear rates in a 
drilling mud 22 1 . 



Above ay, G' drops dramatically and rapidly becomes 
negligible when compared to G" . In other words, the 
sample becomes fluid above ay. 



Following previous works [171 . 122l | , we anticipate that 
the yield stress corresponds to the value of the shear 
stress on the plateau observed during the downward 
sweep. This yields ay ~ 51 Pa in satisfactory agree- 
ment with the previous estimate of 47 Pa, so that one 
has CTy = 49 ± 2 Pa. Finally, the shape of the flow curve 
at low shear rates (7 < 0.5 s~^) points to significant 
wall slip which is highly probable in our smooth geom- 
etry 17|. However, as we shall see below through time- 
resolved measurement of the local velocity, interpreting 
such a non-stationary flow curve is difficult due to com- 
plex temporal behaviours that take place over time scales 
of about one hour. 



Combined velocimetry and rheology 

In this section, we present the combined rheological 
and velocity profile measurements performed under im- 
posed shear rate on a 3 wt. % Laponite sample seeded 
with 1 wt. % microspheres in a smooth Couette cell. As 
recalled in the introduction, our previous experiments 
5J| with a lower microsphere concentration of 0.3 wt. % 
showed that slippery boundary conditions may lead to a 
rather complex scenario of slow fragmentation and ero- 
sion of the initially solid material. In paragraph , we 
shall explore the influence of the imposed shear rate on 
this scenario for a given microsphere concentration of 
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1 wt. % through experiments performed over typically 
5000 s. The influence of the microsphere concentration 
will be addressed in paragraph . 

In paragraph below, we focus on the velocity profiles 
recorded within the first few minutes after the inception 
of shear. We show that the flow may first seem com- 
patible with a simple shear localization picture of critical 
shear rate 7c — 125 s~^. However, as explained later in 
paragraph , shear localization does not correspond to the 
stationary state reached by the system and the experi- 
ments have to be conducted over time scales of the order 
of one hour to probe the actual long-time flow behaviour. 



Early stage after shear start-up 

Figure [3] displays the velocity profiles v{r) averaged 
over t = 30-40 s after various shear rates 7 ranging from 
1 to 500 s~^ are applied at time t = 0. The time win- 
dow for averaging the velocity profiles was chosen so that 
fast initial transients due to flow inception have died out. 
Such transients last typically a few seconds, after which 
the dynamics becomes much slower and the velocity pro- 
files do not significantly change over 10 s. In the fol- 
lowing, we shall analyze these velocity profiles as if they 
corresponded to steady measurements, keeping in mind 
that, over time scales longer than a few minutes, a much 
more complex behaviour will emerge. 

For 7 < 10 s""'^, total wall slip is observed and the 
material undergoes solid-body rotation. In this case, the 
shear rate effectively experienced by the material van- 
ishes, a feature that was already mentioned in Ref. [sst . 
Above 7 = 10 s^^, the velocity profiles are characterized 
by a flowing zone close to the rotor that coexists with 
an "arrested" solid-like region where the local shear rate 
is zero. The size of the fluid-like region increases with 
the imposed shear rate. For 7 > 125 s~^, velocity pro- 
files are homogeneous and the whole sample is fluid on 
the time window investigated here, i.e. after a few 10 s. 
These measurements are reminiscent of shear localization 
as reported in various other SGM [S [22, 23, 26,] . 

Note that apparent wall slip remains significant at both 
walls as long as some "arrested" region is present. In- 
deed, for 7 < 125 s~^, the velocity of the sample in the 
close vicinity of the cell boundaries never reaches that of 
the walls, i.e. v{r = 0) < vq at the rotor (where vq is 
given by Eq. ([T])) and v{r = e) > at the stator. Wall 
slip becomes negligible only when a homogeneous flow is 
recovered for 7 > 125 s~^. In the following, in order to 
compare our data with predictions and experiments ob- 
tained in the absence of wall slip, we shall focus on the 
effective global shear rate 7off defined by [5^: 



7off = 



Rl 



v{Q)~v{e) 



Ri{Ri+R2) e 
rather than on the applied shear rate 7. 



(2) 



By using quadratic fits of the velocity profiles within 
the fluid-like band, one may easily extract the width Vc 
of the flowing zone as well as the local shear rate 710c 
averaged over this sheared region. Such an analysis is 
presented as a function of the effective global shear rate 
7eft in Fig.mja) and (b). The error bars in Fig.Hfb) corre- 
spond to the standard deviation of the fitted local shear 
rate. It can be seen that the proportion of the sheared 
region increases roughly linearly with 7off, consistently 
with the prediction of the simplest theoretical scenario 
for shear localization [lol . [lij and with recent experimen- 
tal findings on a colloidal gel [2^. Here, Fig.|4Ka) clearly 
points to a critical shear rate 7^ ~ 125 s"^ . 

Yet, Fig. Hl^b) contradicts such a simple shear local- 
ization scenario. Indeed, it is clearly seen that the local 
shear rate 710c is not constant and equal to 7c. It rather 
increases sharply from 710c — 60 s~^ to 125 s~^ as 7 is 
increased. This is confirmed by looking at the "local" 
flow curves displayed in Fig. [3]Jc) and showing the local 
shear stress cr{r) plotted against the local shear rate 7(r). 
(j{r) is computed from the global shear stress ct measured 
by the rheometer simultaneously to the velocity profiles 
using [5^: 



a{r) = 



o p2 7j2 



(3) 



where the proportionality factor accounts for the Cou- 
ette geometry. The local shear rate 7(r) = (f ) 

is simply estimated from the derivative of the quadratic 
fits of the velocity profiles of Fig. [31 The resulting local 
flow curves cr(r) vs 7(r) are compared to the downward 
sweep of Fig. in Fig. For the highest shear 

rates, the local data collapse rather well on the global 
flow curve. This indicates that for 7 > 200 s~^, the ma- 
terial becomes fully fluid-like almost instantly and that 
no time-dependent phenomena further occur. For smaller 
shear rates, the local flow curves reveal two interesting 
features: (i) the existence of a "forbidden" range of local 
shear rates since 7(r) < 60 s~^ is never observed and 
(ii) the relevance of time-dependent, history-dependent, 
or metastable phenomena since the local shear stress for 
t = 30-40 s is significantly smaller than the shear stress 
recorded during the downward sweep of Fig. Hfc). Fi- 
nally, for shear-localized velocity profiles, 7(r) covers the 
range 60-125 s^^. 

In the next paragraph, we shall see that shear local- 
ization as revealed in Fig. [3] is only transient so that the 
non-standard scenario shown in Fig. Sfb) and (c) is not 
so surprising. In any case, it is important to stress the 
fact that, if the velocity measurements had been stopped 
after a few minutes, shear localization could have been 
wrongly interpreted as the steady-state for yielding in 
this SGM. 
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Long-time flow behaviour 

As shown in Fig. [SJ the long-time behaviour of our 
sheared samples is not shear localization. A complex 
spatio-temporal scenario rather develops on time scales 
ranging from a few minutes to several hours depending 
on the applied shear rate. As in Ref. [54|, we propose to 



distinguish three different regimes in the temporal evo- 
lution of the flow. 

In a first step called regime I (see Fig. [5ja) and (b)), 
velocity profiles are seen to evolve from shear localization 
towards plug-like flow with strong wall slip within about 
5 min for 7 = 110 s^^. Indeed, the arrested region close 
to the stator "detaches" from the fixed wall and expands 
across the gap, which leads to a velocity profile that is 
flat (and referred to as "plug-like" in the following) over 
most of the cell except for two small sheared regions of 
width ~ 0.15 mm at both walls. At the same time, di- 
rect observations of the sample indicate that the sample 
rotates as a solid body and presents large heterogeneities 
of typical size 1 cm with fracture-like streaks. At the 
end of regime I, the velocity profiles are symmetrical and 
the solid material in the center rotates at roughly wo/2, 
where vq is the rotor velocity. 

In regime II (see Fig. WLc) and (d)), the velocity pro- 
files oscillate within a few seconds between the plug-like 
profile described above and a homogeneous Newtonian- 
like velocity profile with much smaller wall slip (that we 
shall call a "linear" velocity profile below) . This is clearly 
evidenced on the spatiotemporal diagram of v{r,t) in 
Fig. [6l^B2) and on the velocity signals v{rf),t) shown in 
Fig. Eljb) for two different positions in the gap. More- 
over, movies of the sheared samples show that the solid 
material breaks up into smaller pieces surrounded by flu- 
idized zones in both the azimuthal and vertical directions 
[5^ . Simultaneous observations and velocity measure- 
ments also reveal that linear profiles are recorded when- 
ever a homogeneous fluidized region passes in front of the 
USV transducer whereas plug-like profiles correspond to 
solid pieces fioating in the sample [54| . 

Therefore, the temporal oscillations of the velocity 
profiles are simply a consequence of the breaking-up of 
the sample into a very heterogeneous pattern of smaller 
blocks that rotate as solid bodies within the fluidized 
material. Since the time required for recording a single 
velocity profile is much smaller than the rotation period 
of the solid blocks, the proportion $(t) of plug-hke fiow 
profiles recorded during a given time window centered 
around time t can be directly interpreted as the propor- 
tion of solid pieces within the sample. Yet, this obviously 
assumes that the time evolution is the same all along the 
vertical direction, which is not truly the case. In partic- 
ular, it is often observed that the top and bottom of the 
Couette cell get fluidized more rapidly than the center of 
the cell where our velocity measurements are performed. 



This can most probably be attributed to end effects and, 
in the absence of any two-dimensional information on the 
velocity field, we shall neglect the dependence of the fiow 
along the vertical direction. In practice, we discriminate 
between "plug-like" and "linear" profiles by monitoring 
the local shear rate in the middle of the gap: a velocity 
profile is called plug-like when 7(r = 0.5) < 7/2. ^{t) 
is then defined as the fraction of plug-like profiles over a 
time window centered at t and of duration ranging from 
10 to 50 s depending on the applied shear rate. 

In most cases, $(t) decreases slowly throughout regime 
II, which we interpret as an erosion of the solid pieces by 
the fiuidized material through viscous stresses. This is 
directly confirmed by optical imaging that shows a simul- 
taneous decrease of the average size of the solid blocks 
jsot , although we were not able to extract more quantita- 
tive information about the evolution of solid content due 
to the low contrast of the images. In the case of Fig. [5l 
the erosion process takes about 400 s. Figure [6ljB) gath- 
ers the full data set for 7 = 110 s~^ by showing the shear 
stress response cr(t) recorded by the rhcometer, a spa- 
tiotemporal diagram of the velocity profiles v{r,t), and 
$(t). In particular, Fig. [6l^B2) reveals that the temporal 
variations of the local velocity get slower and slower as 
$(i) decreases and progressively give way to a homoge- 
neous stationary flow, which corroborates our picture of 
solid blocks that get smaller and smaller due to erosion. 
Note that such erosion mostly occurs in the azimuthal 
and vertical directions since velocity profiles indicate that 
the radial size of the solid pieces remains comparable to 
the gap width all along regime II. 

Regime II ends when the sample is fully fluidized and 
only linear velocity profiles are recorded (<& = 0). For 
longer times, the sample remains fiuid and no shear lo- 
calization is observed in the steady state called regime 
III (see Fig.jSjc) and (f)). Slip velocities are always very 
small in regime III. 

In summary, regime I corresponds to a transition from 
shear localization to plug-like fiow at vo/2. In regime 
II, the solid region that occupies a large part of the gap 
breaks up into solid blocks that get slowly eroded. In 
regime III, the sample has finally reached its steady-state 
characterized by almost linear velocity profiles with neg- 
ligible wall slip. 



Influence of the applied shear rate 

An important issue concerns the robustness of the frag- 
mentation and erosion scenario described above. Here, 
we explore the infiuence of the applied shear rate. Fig- 
ure [6] shows the results obtained for 7 = 100, 110, and 
120 s^^. In all cases, regimes I, II, and III can always 
be found and are indicated using dashed lines. The most 
striking effect of a rather small increase of 7 is the very 
strong speeding up of the whole process. While it takes 
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about an hour to fully fluidize the sample for 7 = 100 s~^, 
an almost homogeneous flow is recovered within about 
2 min for 7 = 120 s~^. We shall discuss this feature in 
more details in Section . 

Moreover, the stress response a{t) may be analyzed in 
view of the local measurements. Indeed, regime I seems 
to be associated with rather slow stress variations with 
large amplitudes of about 2 Pa. This is most clearly seen 
in Fig.[B{Bl) and (CI), although no systematic behaviour 
of a(t) can be extracted. On the other hand, in regime 
II, a{t) always fluctuates much more rapidly and with 
a much smaller amplitude around an average value that 
slowly decreases. We presume that these noisy features 
in a(t) are the signature of the local fragmentation and 
erosion process. However, due to the highly heteroge- 
neous nature of the flow in regime II, most of the local 
features average out in the global stress response, which 
therefore presents only small fluctuations around a mean 
value close to the yield stress ay = 49 ± 2 Pa. This is 
shown in more detail in Fig. [7] where no clear correlation 
between the stress response and the local velocity can 
be detected. Finally, when the steady state is reached 
in regime III, the shear stress presents a constant value 
with negligible fluctuations. 



Influence of the microsphere concentration 

The fact that our samples contain glass microspheres 
in a rather large amount raises legitimate questions about 
the influence of such seeding particles on our results. To 
investigate the influence of the microsphere concentra- 
tion, experiments were repeated with a smaller amount 
of hollow glass spheres of 0.3 wt. %. Figures [8l [H andflOl 
show that the flow behaviour unveiled above is also ob- 
served for 0.3 wt. % microspheres, as reported in Ref. [s^ l 
where only one applied shear rate (7 = 17 s~^) was in- 
vestigated. In particular. Figs. [S] and [5] are strikingly 
similar. The only differences lie (i) in the lower contrast 
of the images when only 0.3 wt. % microspheres are used 
and (ii) slightly larger slip velocities in the case of the 
lower microsphere concentration. 

The comparison between Figs. [6] and [9] also proves 
the robustness of the fragmentation and erosion scenario. 
Stress responses display broadly the same features (ex- 
cept for the highest applied shear rate in Fig. [9] for which 
an increasing a{t) was recorded). Both the qualitative 
evolutions of v{r,t) and of <I>(t) do not seem to depend 
on the microsphere concentration. 

However, the values of the shear rate at which this 
complex yielding behaviour is observed are noticeably 
different. Indeed, for 1 wt. % microspheres, the frag- 
mentation and erosion process is observed for 7 ~ 85- 
250 whereas this range is 7 ~ 15-65 s^^ for 0.3 wt. % 
microspheres. This remarkable effect of the microsphere 
concentration is further discussed in the next Section. 



Moreover, as revealed by Fig. [10] for a Laponite suspen- 
sion seeded with 0.3 wt. % microspheres and sheared at 
7 = 29 s~^, oscillations may be detected in the stress 
response with the same period as the velocity oscilla- 
tions. Although the relative stress fluctuations remain 
small, such a correlation hints at a less heterogeneous 
flow pattern and possibly more collective fragmentation 
and erosion processes. 

DISCUSSION AND MODEL 
Comparison with other works 

The data presented in this paper show that the frag- 
mentation and erosion scenario first reported in Ref. |54l | 
is found over a large range of shear rates and is robust 
against a change of microsphere concentration. To our 
knowledge, only two other works on SGM close to yield- 
ing may have reported behaviours similar to the one pre- 
sented here. First, intermittent jammed states were ob- 
served together with the occurence of strong wall slip in 
aglassy colloidal star polymer using NMR velocimetry 
2ll |. Second, more recently, oscillating velocities were 
reported in 3 wt. % Laponite samples under low shear 
rates [1^. In the latter case, the two extreme velocity 
profiles were reconstructed thr oug h pointwise DLS mea- 
surements. Figure 4 in Ref. [53| is broadly consistent 
with our observations, e.g., with our Figs.[5Kc) andlDJc). 
Although the oscillating velocity profiles were interpreted 
in terms of a stick-slip instability by the authors, we be- 
lieve that this behaviour may actually be due to fragmen- 
tation and fluidization as described here. Unfortunately, 
in Ref. (Hsj . the evolution of the sample on very long 
time scales was not investigated and optical imaging was 
not possible since the sample did not contain any seeding 
particles and was thus perfectly transparent [60|. 



Characteristic time scale for erosion 

To further discuss these results, we extracted a char- 
acteristic time T for the erosion process defined as the 
time it takes for $(t) to relax from 0.9 to 0.1 during 
regime II (see also the solid lines in Figs. [n]and[5]). In 
Fig. fTlT a'). this characteristic time is plotted as a function 
of the applied shear rate for the two microsphere concen- 
trations investigated here. It clearly appears that r varies 
by more than three orders of magnitude within a rather 
small range of shear rates of order 20 s~^. The time scale 
for erosion seems to diverge for shear rates smaller than 
7* ~ 85 s""'^ for 1 wt. % microspheres and 7* ~ 15 for 
0.3 wt. % microspheres. In fact, for 7 < 7*, only plug-like 
flow profiles were observed over the longest experimen- 
tal waiting times of about 10^ s, so that fragmentation 
actually never occured. For the largest shear rates, the 
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material becomes fluid within seconds and, due to our 
limited time resolution of about 1 s per velocity profile, 
relaxation times shorter than a few seconds cannot be 
accessed. 

The large scatter of these data may be attributed to 
poor reproducibility of the yielding behaviour. Indeed, 
experiments performed on two different samples from the 
same batch at the same shear rate may yield estimates for 
T that differ by almost a decade. Uncertainty also arises 
from the difficulty in precisely locating the beginning and 
the end of regime II. As seen for instance in Figs. ^A) 
and[9|B), ^{t) may start to decrease and go back to $ = 
1 within a few minutes so that fragmentation seems to 
begins but soon stalls before starting again. Figure [6lJC) 
also shows that very small solid pieces may persist for 
very long times even though <I>(t) drops to values of order 
0.2 in less than 10 s. Such behaviours, which could also 
be ascribed to spatiotemporal heterogeneities of the flow 
in the vertical direction, lead to large error bars on t. 

Still, the very steep decrease of r for 7 > 7* prompts us 
to model the experimental data by the phenomenological 
law: 



A(7-7*)" 



(4) 



where n is some characteristic positive exponent. Due 
to the large spread of the data, fits to equation ^ with 
three free parameters do not yield reliable results. How- 
ever, forcing the value 7* = 85 s~^ and excluding the 
most deviant data points leads to a satisfactory fit with 
n = 2.1 and A = 2.2 10^ for 1 wt. % microspheres (see 
solid line in Fig. [TlT a)). Note that when 7* is varied 
from 75 to 95 s~^, the exponent n varies from 2.5 to 1.5. 
Therefore, in order to account for the large experimen- 
tal scatter, we shall conclude that 7* = 85 ± 10 s"-'^ and 
n = 2.1 ± 0.5 provide a good modelling of the 1 wt. % 
microsphere data. 

In Fig. ITlT b). T is plotted in logarithmic scales against 
7 — 7* using 7"^ = 85 s~^ for 1 wt. % microspheres and 
7* ~ 14 s~^ for 0.3 wt. % microspheres respectively. The 
rather good collapse of both data sets allows us to con- 
clude that n and A do not significantly depend on the 
microsphere concentration while the two values for 7* 
are clearly different. 



A simple model for erosion 

In order to model the erosion process, we suppose that 
the sample at the beginning of regime II is characterized 
by some distribution of solid fragments of initial size Rq. 
To account for the progressive fragmentation of the sam- 
ple at the end of regime I, each fragment is assumed to 
start being eroded at an initial time to characterized by 
a given probability distribution. 

Erosion results from viscous friction at the surface of 
the solid pieces. Since the solid pieces have a yield stress 



tTy, the viscous stress available for erosion is 77/7 — ay, 
where r]f is the viscosity of the surrounding fiuid, here, 
the fluidized Laponite suspension. Rewriting this stress 
as 77/(7 — 7"*), with 7"* = ay/rjf, shows that the charac- 
teristic time scale for erosion is 1/(7 — 7"^). 

We then introduce the size rt of the "elementary brick" 
that detaches from solid fragments due to erosion. Since 
erosion occurs at the surface, we expect the number of 
bricks = R{t)^/r^ in a fragment of size R{t) at time 
t to be linked to the number of bricks at the surface 
= Rit)yrl by: 



dt 



oc -(7 - 7*)A^s 



This simple dimensional argument leads to: 



dR 
'dt 



(5) 



(6) 



up to some dimensionless multiplicative constant of order 
unity which we may incorporate into rt,. Note that the 
same equation for R(t) would be found if one rather con- 
siders a two-dimensional problem where the solid blocks 
are discs of constant thickness, which may be more rele- 
vant for the present experiments. 

Finally, we compute ^{t) from the distribution of frag- 
ment sizes {^j(i)}j=i...jv at time t as: 



(7) 



where is the total number of fragments and Rjq is the 
initial size of fragment number j whose erosion starts at 
time tjo. From Eq. one has: 



Rjo if t < tjo 
for i > tjo 4 



(8) 



R 



■jo 



''6(7 ~ 7*)(^ ^ tjo) otherwise. 



is the total erosion time for frag- 



where Ti = . ^. ^ 

ment number j. 

The definition of <I>(t) in Eq. ([7]) was chosen in order 
to be consistent with the experimental measurements of 
$(i) presented in Section . Indeed, these measurements, 
based on the proportion of plug-like velocity profiles mea- 
sured within a given time window, correspond to one- 
dimensional measurements and should be related to the 
linear size Rj{t) of the solid pieces. Furthermore, al- 
though the exact shape of $(t) depends on the details of 
the distributions chosen for Rq and to, it is clear from 
Eq. ([5]) that the characteristic time scale for ^{t) should 
be of the order of: 



n 7-7* 

where Ro is the average initial fragment size. 



(9) 
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Comparing with Eq. ([4]), our simple model predicts 
a n = 1 scaling instead of the experimental exponent 
n = 2.1 ± 0.5. This is most probably due to the crude as- 
sumptions (i) that the distributions for Rq and Iq (in par- 
ticular i?o) do not depend on the applied shear rate, (ii) 
that erosion is governed by the detachment of elementary 
bricks whose size is independent of both time and shear 
rate, and (iii) that there is a direct link between the yield 
stress of the bulk material at rest and the critical shear 
rate 7*. Still, the model predicts a rather good order of 
magnitude for r. Indeed, if r?, is estimated from the elas- 
tic modulus Go as n ^ (fcsT/Go)^/^ with Go = 500 Pa, 
one finds rf, ~ 40 nm, which roughly corresponds to the 
diameter of the Laponite discs. Taking i?o '--^ 1 cm for 
the initial size of the fragments as observed in the exper- 
iments and 7 — = 10 s~^, Eq.® leads to r ^ 10'* s in 
statisfactory agreement with the measurements reported 
in Fig. [m 

In order to illustrate our simple model, we chose to 
focus on the case where Rjo is fixed to a given value 
i?o for all solid fragments and tjo is characterized by a 
Gaussian distribution of mean to a-nd standard deviation 
dto/V2. In this case, for TV > 1, Eqs. ^ and dH) become 
respectively: 



Rq 



p{to) dto , 



(10) 



where p{to) = exp(— (to — to)"^ / Stg) / (y^Sto) is the prob- 
ability density function of the initial time to, and 



Rf) 



if to < t - r 
if t - T < to < t 
if to > t , 



(11) 



where we have noted r 



Ro 



the erosion time for 



rt T-7* 

a single fragment. Inserting Eq. (fTTj) into Eq. pO|) and 
using the dimensionless time t = (t — to)/(5to leads to the 
following analytical expression: 



$(t) = 



1 — erf(t — a) 



1 - 



t-y 



(12) 

where we have defined the ratio of the erosion time to the 
standard deviation of to as a = r/Sto and "erf" stands 
for the error function: crf(x) = e~^' dy. 

Figure [T2l shows the predicted $(t) for various values 
of a. In the limit a 1, where the distribution of initial 
times is very wide compared to the erosion time, Eq. 
reduces to $(t) = (1 — erf(t))/2 as shown by the red 
dashed line. On the other hand, the opposite limit a 3> 1 
corresponds to very slow erosion so that the initial times 
can all be considered as equal to to and the dynamics of 
$(t) is simply that of a single fragment, i.e. a linearly 
decreasing function (see black dashed line). 

Finally, a comparison between the predictions of 
Eq. (|12p and experimental $(t) data taken from Figs.lH] 



and [9] is shown in Fig. [131 In the model, the average 
time ti is directly estimated from the experimental data 
by looking for the time at which regime II begins. For a 
given microsphere concentration, good agreement could 
be found for two different applied shear rates by using 
the same standard deviation Sto and varying only the 
erosion time r. The transition from linear to curved 
shapes for $(t) as the erosion time scale is decreased, 
i.e. as the shear rate is increased, constitutes the most 
prominent feature of our simple approach. In spite of 
large experimental scatter, this feature is noticeable at 
least in Figs. fTSFa) and (d). However, as already noted, 
the model is not compatible with the strong divergence 
of r as 1/(7 — 7*)^. Within the erosion model proposed 
above, the experimental scaling thus points to an addi- 
tional dependence of the initial fragment size Rq on cis 
1/(7-7*)- 



Open questions 

The main open issue raised by our results concerns the 
mechanism for fragmentation and its relation to bound- 
ary conditions. As reported in Ref. [H^l, fragmentation is 
only observed with "smooth" surfaces and must be some- 
how linked to boundary conditions. We proposed to ex- 
plain the observed fragmentation by temporary sticking 
events localized at specific sites where the yield stress 
may be overcome, giving rise to a fluidized region sepa- 



rating solid pieces [54| . Note that such events are highly 



probable because the wall roughness (about 15 nm for our 
"smooth" cell) is comparable to the size of the Laponite 
platelets. 

However, the fact that the characteristic shear rate 7* 
strongly depends on the amount of seeding microspheres 
calls for a more subtle interpretation. One may even 
wonder whether the fragmentation and erosion scenario 
would still be observed if the microsphere concentration 
was further reduced or if the Laponite suspension did not 
contain any microsphere at all. The observation of oscil- 
latory behaviours in pure Laponite 5^, G^] tends to sup- 
port the possibility of our scenario even in the absence of 
any seeding particle. In our experiments, we lowered the 
microsphere concentration down to 0.1 wt. %, which still 
allowed us to measure velocity profiles with reasonable 
accuracy. Fragmentation and erosion could be observed 
in our "smooth" cell only when wall slip was promoted 
in an even better way by lubricating the cell walls using 
vacuum grease (data not shown). We also tried to seed 
our Laponite suspensions with polystyrene (rather than 
glass) microspheres. This led to the same observation 
that the fragmentation and erosion scenario needed very 
slippery boundary conditions to develop. 

Therefore, we may speculate that surface roughness 
is not the only relevant control parameter but that 
microspheres-Laponite and/or microspheres-surfacc in- 
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teractions may also play a significant role and influence 
the values of 7* and r. However, a systematic study of 
these effects, as well as those of sample age and Laponite 
concentration, still remains to be performed. 

Finally, other important open questions concern (i) the 
link between the critical shear rate 7* for erosion and 
the characteristic shear rate 7c involved in shear local- 
ization at early times or in the absence of wall slip and 
(ii) the possible effects of confinement on the observed 
shear-induced fragmentation. 

CONCLUSION 

In this work, wc have investigated in details the origi- 
nal yielding scenario of Laponite suspensions seeded with 
glass microspheres first reported in Ref. (H^j and trig- 
gered by slippery boundary conditions. Such a scenario 
involves shear localization at short times followed by a 
slow evolution towards plug- like flow (regime I). Then, 
over much longer time scales, fragmentation and erosion 
are observed provided the applied shear rate is larger 
than some critical value 7*. Fragmentation leads to ap- 
parent oscillations in the velocity profiles (regime II). Ex- 
cept for small noisy and/or oscillatory features in the 
stress response, the very large local fluctuations of the 
flow field have no significant impact on the global rheo- 
logical response. Depending on the shear rate, a homoge- 
neous flow (regime III) is recovered within a few minutes 
up to several hours. This scenario was reported for a 
wide range of shear rates and for two different micro- 
sphere concentrations. The characteristic time scale for 
erosion roughly follows a 1/(7 — 7*)^ power law, where 
7"^ depends on the microsphere concentration. A sim- 
ple model was introduced to qualitatively reproduce the 
erosion process. Yet, a better understanding and charac- 
terization of the fragmentation and erosion phenomena 
in such a soft glassy material are still needed both from 
the experimental and from the theoretical sides. These 
results also prompt us to look for similar yielding be- 
haviours in other materials under slippery conditions and 
on long time scales. 
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FIG. 1: Sketch of the experimental setup for combined rhe- 
ological and velocity profile measurements. The Couette cell 
is shown in gray. 




FIG. 2: Rheology of a 3 wt. % Laponite suspension seeded 
with 1 wt. % microspheres. Viscoelastic moduli G' (□) and 
G" (•) as a function of (a) oscillation frequency for a fixed 
shear stress amplitude of 0.5 Pa and (b) shear stress amplitude 
for a fixed frequency of 1 Hz. (c) Flow curve cr vs -y obtained 
by sweeping up (>) the shear rate for 200 s and down (<l) for 
200 s. All measurements are performed on fresh samples from 
the same batch using the protocol described in Section . 
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FIG. 3: Velocity profiles recorded during the early stage of 
our experiments on a 3 wt. % Laponite suspension seeded 
with 1 wt. % microspheres for imposed shear rates 7 — (a) 
1 (green), 5 (blue), 11 (red), 20 (black), (b) 30 (green), 40 
(blue), 50 (red), 60 (black), (c) 70 (green), 80 (blue), 100 
(red), 125 (black), (d) 150 (green), 200 (blue), 300 (red), and 
500 s"'^ (black), r is the radial distance from the rotor. The 
velocity profiles were averaged over t — 30-40 s. The solid 
lines show quadratic fits of the data performed over the flow- 
ing region. 
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FIG. 4: (a) Width and (b) local shear rate 7ioc of the 
fluid-like region as a function of the effective global shear rate 
7off- (c) Local rheological flow curves a{r) vs 7(r) (colors) 
together with the global rheological data recorded during the 
downward sweep shown in Fig. [2](<). The colors correspond 
to the various velocity proflles of Fig. |3] (see text). The dotted 
lines in (a) and (b) indicate shear localization with a constant 
> = 125 s-^ 
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FIG. 5: Long-time behaviour of a 3 wt. % Laponite suspen- 
sion seeded with 1 wt. % microspheres for 7 = 110 s~^. Ve- 
locity profiles taken at times (a) t = 1 (■), 185 (o), 222 {<), 
and 302 s (>) in regime I, (b) t = 478 (A) and 493 s (t) in 
regime II, and (c) t = 1005 s (♦) in regime III. The velocity 
profiles v{r) have been rescaled by the rotor velocity vq. The 
right column shows pictures of the sample typical of the three 
regimes. The height of each picture is 30 mm. 
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FIG. 6: (1) Shear stress response (T(t), (2) velocity profiles 
v{r,t), and (3) fraction of plug-like velocity profiles in 
a 3 wt. % Laponite suspension seeded with 1 wt. % micro- 
spheres for (A) 7 = 100 s"\ (B) 7 = 110 s~\ and (C) 
7 = 120 s~^. The spatiotemporal diagrams in (2) use lin- 
ear gray scales ranging from v = (black) to « = 7e (white) . 
The solid lines in (3) are exponential relaxations shown to 
guide the eye across regime II. 
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FIG. 7: Zoom into regime II for a 3 wt. % Laponite suspension 
seeded with 1 wt. % microspheres under 7 — 100 s~^. (a) 
Shear stress response a{t) and (b) local velocity v{ro,t) for 
ro =0.1 mm (blue) and ro — 0.8 mm (red). 
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FIG. 8: Long-time behaviour of a 3 wt. % Laponite suspension 
seeded with 0.3 wt. % microspheres for 7 = 39 s^^. Velocity 
profiles taken at times (a) t = 14 (■), 178 (o), 955 and 
1399 s (>) in regime I, (b) t = 1933 (a) and 1974 s (▼) in 
regime II, and (c) t — 2857 s (♦) in regime III. The velocity 
profiles v{r) have been rescaled by the rotor velocity vq. The 
right column shows pictures of the sample typical of the three 
regimes. The height of each picture is 30 mm. 
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FIG. 9: (1) Shear stress response (T(t), (2) velocity profiles 
v{r,t), and (3) fraction of plug-like velocity profiles in 
a 3 wt. % Laponite suspension seeded with 0.3 wt. % mi- 
crospheres for (A) 7 = 29 s"\ (B) 7 = 39 s"\ and (C) 
7 = 50 s~^. The spatiotemporal diagrams in (2) use linear 
gray scales ranging from v = (black) to v = -ye (white). The 
solid lines in (3) are exponential relaxations shown to guide 
the eye across regime II. 
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FIG. 10: Zoom into regime II for a 3 wt. % Laponite suspen- 
sion seeded with 0.3 wt. % microspheres under 7 = 29 s~^. 
(a) Shear stress response a{t) and (b) local velocity v{ro,t) 
for To = 0.1 mm (blue) and ro = 0.8 mm (red). 




FIG. 11: Characteristic erosion time r in regime II vs (a) the 
applied shear rate 7 and (b) 7 — 7* for 3 wt. % Laponite 
suspensions seeded with 1 wt. % (□) and 0.3 wt. % (•) mi- 
crospheres. The solid lines are r = ^(7 — 7*)^" with n = 2.1, 
A = 2.2 10^ and 7* = 85 s"^ (7* = 14 s"^ resp.) for the 
1 wt. % (0.3 wt. % resp.) microsphere concentration. 




(t-V/6t„ 

FIG. 12: Fraction of solid fragments <I> computed from 
Eq. (|12|l as function of the reduced time t — {t — to)/5to 
for different values of a — r/Sto (see text). The red and 
black dashed lines are respectively $(i) = (1 — erf(f))/2 and 
$(?) = 1 - r/a for < r< a = 10. 
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FIG. 13: Predictions of Eq. (|12|l (thick solid lines) compared 
to experimental data (o). (a) to = 210 s, Sto — 20 s, and 
T = 360 s (a = 18). (b) = 110 s, Sto = 20 s, and t = 10 s 
{a = 0.5). (c) IE = 1430 s, Sto = 200 s, and r = 1600 s 
{a = 8). (d) U = 1750 s, Sto = 200 s, and r = 400 s (a = 2). 
The experimental conditions, namely the microsphere weight 
fraction and the applied shear rate, are indicated directly on 
the figures. 



